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Abstract 

We study solutions of the decoupled Maxwell equations in the exterior region of a Schwarzschild black 
hole. In stationary regions, where the Schwarzschild coordinate r ranges over 2M < n < r < r2, we 
obtain a decay rate of for all components of the Maxwell field. We use vector field methods and do 
not require a spherical harmonic decomposition. 

In outgoing regions, where the Regge- Wheeler tortoise coordinate is large, r* > et, we obtain decay 
for the null components with rates of \4>+\ \a\ < Cr~^^'^, \(j)o\ ~ |p| + |o"| < Cr~^\t — rf\~^^^, and 
|0_i| ~ |q| < Cr~^\t — Along the event horizon and in ingoing regions, where r* < 0, and when 

i + r, > 1, all components (normalized with respect to an ingoing null basis) decay at a rate of Cu+~^ 
with 11+ = t + r» in the exterior region. 

1 Introduction 

The subject of this paper is the study of decay of solutions to the decoupled Maxwell equations in the exterior 
of a Schwarzschild black hole. The Maxwell field is a 2-form which we may write in abstract index notation 
as an antisymmetric (0, 2)-tcnsor field on a manifold A4, 

Fen^iM) or F^p^-Fp^. 

It satisfies the Maxwell equations: 

*d*F=0 or \7°'Fafj=0 (1) 

dF =0 V[„i^^^] =0. (2) 

The exterior region of the Schwarzschild solution is a Lorentz manifold on which the metric is given in terms 
of coordinates t G M, r > 2M, (9, (j)) e 5^ by 

ds^ = - (1 - 2M/r)dt^ + (1 - 2M/r)-^dr^ + r^{d0^ + sin^ Odcj)^). (3) 

This problem comes from general relativity. In general relativity, a model of the universe consists of a 
space-time manifold A4, possibly fields describing matter, and a Lorentz (pseudo-) metric g which satisfies 
Einstein's equation. Gravity is described by the curvature of g. The simplest and longest-known solution is 
Minkowski space, R^+'^ with the flat metric —dt^ + dx^ + dy'^ + dz'^ . After this, the Schwarzschild manifold is 
the longest-known solution to Einstein's equation. It is the paradigmatic example of the class of black hole 
solutions, which play an important role in relativity. The Maxwell field describes electromagnetic radiation. 



In Einstein's equations, the energy-momentum tensor of the matter fields should influence the curvature. 
By decoupled, we mean that the electromagnetic field does not influence the Schwarzschild solution, which 
is taken as a fixed background manifold. We call the Schwarzschild solution the Schwarzschild manifold and 
use the word solution to refer to solutions of the Maxwell equations 

Since is a tensor, there is no coordinate independent norm with which to measure it (or, at least, not 
all components of it). To discuss the decay of F, we make a choice of basis and show that the corresponding 
components decay. A simple choice of basis consists of the coordinate vector fields rescalcd so that they have 
unit length X)| = 1). The rescaled vectors are 

f ={l-2M/r)-'^/^dt, R={l-2M/ry/^dr, Q =r-^de, i> ^r'^ siniOy^d^,. 

Given a time-like vector, there is a natural decomposition of the Maxwell field into electric and magnetic 
components. Since the Schwarzschild manifold has a time-translation symmetry, this provides a natural 
choice of time-like direction, T. The corresponding electric and magnetic components are 

X e {i?,e,$}, 

X,Y, Z a. cyclic permutation of i?, G), $, 



Now that we have a choice of components for the Maxwell field, it is possible to state the main decay 
result of this paper. 

Theorem 1 (Decay in stationary regions). Let 2M < ri < r2 < oo. There is a constant C and a norr^ 
H[F]{Q) depending only on F and its derivatives on the hyper-surface {0} x (2M, oo) x S"^ such that if F is 
a solution to the Maxwell equations then for all t £ M, r £ [ri, r2], {9, 0) £ S*^, 

\E\ + \B\ <C{l + \t\)-^H[F]{{)). 

The major advance of this work is to find decay rates which govern all components of the (decoupled) 
Maxwell field explicitly. The rates we obtain for stationary regions with r £ (ri, r2) are significantly slower 
than the rate of t~^/'^ which can be obtained in Minkowski space using vector field methods and the rate 
of which was derived formally for the Schwarzschild manifold [22l [23] . Outside of outgoing light-cones, ie 
where t < r» = r -|- ln((r — 2M)/2M) -I- C, decay rates at the same rate as in Minkowski space have already 
been obtained [18j . In the outgoing region, we obtain similar results, which we explain below. Certain 
components of the Maxwell tensor satisfy a scalar wave equation. These components are the zero-weight 
(spinor or null) components. Previous results for wave equations were sufficiently strong to prove decay for 
the zero- weight component with a rate of t~^ in stationary regions and the appropriate decay in outgoing 
regions [1], although this application was not explicitly stated. L"^^ decay without a rate has also been 
explicitly obtained using very different techniques [14| . The existence and asymptotic completeness of wave 
operators taking data on the initial surface t = to the surfaces at r = 2M and at infinity has also been 
shown [1]. 

Our method starts by using the energy-momentum tensor to generate a positive, conserved energy from 
the time translation symmetry and a stronger "conformal energy" from a vector field K. This follows ideas 
in [S] and is very closely related to the analysis of the wave equation in [51 [H [5] . Before the wave estimates 
were known, a similar method was used [18] . The growth of the conformal energy is bounded by a "trapping 
term" consisting of the Ef^ and Bj^ components localized near the photon sphere, r = 3M. In the geometric 

^The norms used are stated explicitly in section |4] For this norm to be finite, it is sufficient that the initial data and its 
first eight derivatives are bounded and decay like r~(^/^+''' (see remark [9)l. The initial data does not need to decay at the 
bifurcation sphere, r — > 2M. We do not use a spherical harmonic decomposition in our analysis; however, from the structure of 
the Maxwell equations, spherically symmetric solutions can have no time dependence and cannot decay sufficiently rapidly for 
the norm H[F]{0) to be finite (see appendix 



Bx =Fyz 

I^p=ii?^p + i4|2 + i4|2, 

\B\'^\B.\' + \B^\' + \B^\\ 
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optics limit, electromagnetic radiation follows null geodesies, which can orbit at r = 3M. Energy can decay 
arbitrarily slowly from this region, at least for the wave equation [24]. Thus, it should be expected that 
there is an obstruction to dispersion near this surface. The trapping term can be controlled because the Ej^ 
and Bj^ components each satisfy a scalar wave equation of a type that's been previously studied [1]. This 
wave equation and the terminology "zero-weight components" for and follow from the analysis of 
the Price equations ([^^ - ([25)1 first appearing in We refer to this reduction to scalar wave equations as 
"spin reduction" . The control on the conformal energy allows us to conclude: 

Lemma 2. There is a constant C and a norm H[F]{0) depending only on F and its derivatives on the 
hyper-surface {0} x (2A/, cx)) x S'^ such that if F is a solution to the Maxwell equations (HI)-©, then for any 
2M < ri < r2 < oo and t sufficiently large, on the surfac^ S ~ {t} x [j'l, 7-2] x S"^ , 

J (\E\'^ + iBp) (1 - 2M/rydr^d^uj <Ct^^H[F]{Of . 

From Soolev estimates and integrated decay estimates, like lemma O for the Lie derivative of F, it is 
possible to prove pointwise decay estimates. In Minkowski space, the four coordinate directions generate 
symmetries, so that the Lie derivatives of a Maxwell field also satisfies the Maxwell equations. Although we 
lack a full set of symmetries, we do have 3 from the time-translation and angular-rotation symmetries. To 
control a fourth direction, we use the Maxwell equations to "trade" the derivatives in the directions of the 
three symmetries for a radial derivative. With Lie derivatives in all directions controlled, we conclude that 
theorem [T] holds. 

To further explain our results and those of others, we describe the geometry of the Schwarzschild manifold 
and its importance. This description can be found in most introductory relativity texts (ie [lTl[20|). The 
Lorentz metric is most simply given in terms of coordinates (t,r,9,(l)) by ([3]). As r — > 00, this metric 
approaches the fiat, Minkowski metric written in spherical coordinates, ds^ = —dt^+dr^+r'^{d9'^ + sin'^9d4P'). 
For r > rg > 2M, the Schwarzschild solution describes the space-time of a vacuum outside a star of radius 
ro and mass M . The restriction on r can be relaxed by considering extensions of this manifold. The metric 
is clearly well-defined in the exterior region i e R, r G (2M, cx)), {9,(f)) G S*^, and in the interior region 
t G K, r G (0, 2M), {0,4)) G . In the interior region, since (1 — 2M/r) is negative, r is a time-like 
coordinate, and t is space-like. The maximal analytic extension of any open subset of the Schwarzschild 
solution is illustrated in the conformal diagram in figure [U in which the angular variables are suppressed. 
There are two exterior regions (/ and ///) and two interior regions (// and IV). By an appropriate choice of 
coordinates, each interior can be smoothly joined to each exterior along a null surface r = 2M . The manifold 
is also smooth at the bifurcation sphere where the four regions meet. However, as r ^ 0, the curvature 
polynomial diverges. 

The Schwarzschild manifold is a prototypical solution to Einstein's equation which has inspired many 
key concepts in general relativity. The asymptotic approach of the metric to the flat, Minkowski metric is 
known as asymptotic flatness. In the conformal compactiflcation of each exterior region of the Schwarzschild 
solution, each outgoing geodesic (with r — > 00 as < — > 00) ends on future null infinity 3"*", and each ingoing 
null geodesic (with r — > 00 as i — > —00) starts on past null infinity 3~ . In essence, an asymptotically flat 
manifold is defined to be one with a future null infinity. A black hole is a region of space-time which cannot 
be joined by future-directed, null or time-like curves to 3^, and an event horizon is its boundary. The future, 
interior region of the Schwarzschild solution is a black hole, and the surfaces where r = 2M are the event 
horizons. In the Schwarzschild manifold, the (future) singularity at r — > is separated from by the event 
horizon. The singularity theorems state that under a broad range of conditions, future singularities must 
form [17| . The weak cosmic censorship conjecture asserts that under some genericity condition, which is not 
yet known, future singularities are always separated from 3"*" by an event horizon. There is a three parameter 
family of asymptotically flat, known, exact solutions to Einstein's solutions which represent massive, rotating, 
charged black holes. This is the Kerr-Newman class, and the Schwarzschild solutions are the solutions with 
positive mass and zero angular momentum and charge. These are stationary, in the sense that they have 

more general surface is permitted in the statement of lemma [5] 
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Figure 1: A conformal diagram for the maximal extension of the Schwarzschild manifold (suppressing the 
spherical coordinates). Thin lines represent boundary points at infinity. Thick lines represent the singularity 
at r ^ 0. Dotted lines represent the event horizon. Regions / and /// are exterior regions, and regions II 
and IV are interior regions. The surfaces 3^ represent future and past null infinity. The points represent 
future and past null infinity. The points i'^ represent spatial infinity. 



a time-translation symmetry sufficiently close to null infinity. These solutions also have singularities, but 
if the angular momentum and charge are beneath a critical threshold, then the singularities are separated 
from the asymptotically flat regions, in the sense that a future-directed, time-like curve from a point in 
an exterioir region will either escape to null or space-like infinity or cross the event horizon, but not both. 
The Kerr-Newman solutions are believed to be the only asymptotically flat, stationary solutions. Physicists 
believe that all black holes should approach one of the stationary, Kerr-Newman solutions. It is not yet 
known if a small perturbation of a Cauchy surface for one of the Kerr-Newman solutions will evolve into a 
solution which remains similar to one of the known solutions. This is the question of black hole stability. 

Stability for Minkowski space was a major and difficult result [5]. Einstein's equations are a complicated 
system of nonlinear equations in which the geometry is dynamic. The linearization of Einstein's equations 
about Minkowski space forms a system called the spin 2 field equations. Obtaining decay estimates for the 
spin 2 field was one step in this proof [5] . Decay estimates for the decoupled Maxwell equations were proven 
at the same time. 

The question of stability of the Schwarzschild solution has also inspired the study of linear fields. In 
the linearization of Einstein's equation, certain components are determined by the solution to a simple wave 
equation [25], and the remaining components arc determined by the solution to a more complicated wave 
equation |29j . Using spinors, Price was able to present a more unified presentation for all components of 
several important, physical systems, according to their spin. Any wave equation is said to have spin 0. 
The Dirac system has spin 1/2. The Maxwell field has spin 1. For any solution of Einstein's equations in 
vacuum, the non-vanishing components of the curvature satisfy certain relations from the Bianchi identities 
and Einstein's equation. In Minkowski space, since the curvature is zero, solutions to the linearization of 
Einstein's equation satisfy the same relations, which are called the spin 2 field equations. In the Schwarzschild 
manifold, since the curvature is non- vanishing, in the linearization of Einstein's equations, there are additional 
terms arising from the derivative with respect to the perturbed metric of the original Christoffel terms. Thus, 
Price distinguishes between the spin 2 field equations and the linearization of Einstein's equation. Formal 
arguments suggest a rate of for fields of all integer spin and for the linearization of Einstein's equations 
[22] . A similar, spinorial presentation of these systems has been made for the Kerr-Newman solutions [27] . 
In each case, certain components were found to satisfy scalar wave equations and then acted as forcing terms 
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in the equations governing the remaining components. 

Most of the subsequent analysis has been focused on the decoupled wave equation. The literature is vast, 
and we list only some of the results. Solutions arc known to remain uniformly bounded in time |19j . The 
scattering theory, concerning the map from the initial data to the the limit on J+ and the event horizon, 
has also been studied on the Schwarzschild manifold and on the more general Kerr-Newman solutions 
|15| . On the Schwarzschild manifold, vector- field techniques have been used to obtain decay results in three 
main steps [4l[2j[9]. First, the vector field K is used to introduce a conformal energy, which is not conserved 
because of trapping. Second, a radial vector field is used to prove a local decay estimate to control the 
trapping term. In R^"'"'^, the radial derivative can be used to make a somewhat similar estimate |21j . In 
K^+'^j estimates involving K and the radial vector field are both referred to as Morawetz estimates. In this 
step, a spherical harmonic decomposition was used in the proofs, but this is no longer necessary [5]. Because 
the scalar equation governing the zero-weight components of the Maxwell field has a simple structure, in 
appendix [B1 we are able to modify the earlier method to obtain decay without using a spherical harmonic. 
Third, the conformal energy is used to control norms. In [11 [9], a strong local decay estimate is proven 
and additional angular derivatives are used to obtain a decay rate of and a similar decay rate in 
outgoing regions. In [S], an additional vector field, y, is used to also prove decay estimates along the event 
horizon. These require weighted or norms of the initial data to be bounded. By an norm, we 
mean, roughly speaking, that the /cth derivative of a solution u is square intcgrable. In [2], only weighted 
H^^'^ norms of the initial data are needed, but a weaker local decay estimate was obtained, which led to less 
control on the conformal energy and a decay rate oit~^/^ for a weighted, spatial norm. Using an entirely 
different technique, based on a representation of the propagator, L"^^ decay has been proven for the wave 
equation on subcritical Kerr-Newman solutions [13] . 

The spin 1/2 system is the Dirac model for the electron. On the sub-critical Kerr-Newman solutions, 
scattering results [16] and L"^^ decay [12] have also been proven. 

For the linearization of Einstein's equations about the Kerr-Newman solutions, the equations found in 
[27] were found to have no unstable modes [28]. For the linearization about the Schwarzschild solution, 
the simpler equations in [5S] satisfy an integrated decay estimate [3] and decay [2]. Although the 
application is not explicitly stated, decay at a rate of follows from (J. 

For the full Einstein equations on a black hole background, decay results are known in the spherically 
symmetric case, when the Einstein equations are coupled to a scalar wave equation and the Maxwell field 
jlOj . By Birkhoff's theorem (see [17]): the Schwarzschild manifold is the only spherically symmetric solution 
to Einstein's equation (treating Minkowski space as the special sub-case with M = 0). The decay rate 
obtained for the scalar field is u+~^^^ along the event horizon. The decay rate is known as Price's 

law and has important implications for the strong cosmic censorship conjecture [7]. 

To discuss decay outside regions of fixed r, it is necessary to introduce components with respect to a 
null tetrad, a basis built from null vectors with certain properties. Physicists may know these as spinor 
components |221 126j. and mathematicians, as the null decomposition [S]. We present one tetrad here and 
discuss exactly what we mean by a null tetrad in section [21 

We start by introducing the Regge- Wheeler radial coordinate, r*, defined by 

dv 

— ^{1 - 2A//r), r(0) =3A/. 

The exterior region of the Schwarzschild solution is given by (i, r, , 0, (p) ranging over M x M x S''^. In these 
coordinates, the Lorentz metric becomes 

g = - (1 - 2M/r)dt^ + (1 - 2Mlr)drJ^ + r^{de'^ + sin^ 9d(f^). 

From this form of the metric, it is clear that any multiple of the vectors dt ± dr, are null. To define our 
null tetrad, wc use ca and cb to denote an orthonormal basis of tangents vectors to 5*^ and e to denote the 
antisymmetric, Levi-Civita tensor on S'^. The null tetrad we will use to state our results is 

{dt+dr,), {l-2M/r)-\dt-dr,), r-^CA, r-^BB. 
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The covariant derivative of this null tetrad along ingoing, radial, null geodesies is zero. Thus, having found a 
natural choice of null tetrad on the initial surface i = 0, we have extended it to the entire future of the initial 
surface by parallel transport along null geodesies falling into the black hole. This is useful for considering 
limits as r ^ 2M. Had we extended the basis by parallel transport along outgoing null geodesies, to study 
the problem as r ^ oo, the factor of (1 — 2M/r)~^ would have been on {dt + dr,) instead of {dt — 9r, )■ 
However, since (1 — 2M/r) 1 as r ^ oo, the difference between our choice of null tetrad and the natural 
choice is vanishingly small. Therefore, we use our choice of null tetrad throughout the future, t > 0. 

The decay of these spinor or null components is not simply a decay in time. This is known from the 
behavior in R^+^. In that case, the heuristic is that the bulk of solutions to the Maxwell equations travel 
out along the light-cone f ~ |a;|. In any fixed region, there is decay because the wave leaves the region. As 
the light-cone expands, the average value of the wave intensity drops. Moving with the wave, the intensity 
decays as it is spread over the increasing area of the light-cone. Thus, some of the decay occurs as a result 
of the wave being far from the light cone, and some occurs as a result of the light-cone being very large. 
Similar behavior occurs on the Schwarzschild manifold. The null coordinates 

=t + and u_ =t — 

are used to measure the distance from the light-cone, and r is used to measure the size of the light-cone in 
the outgoing direction. In the ingoing direction, the radius of the surface of the light-cone also goes like r, 
but since this approaches 2M, the decay occurs only in the null coordinates. 

Theorem 3 (Decay outside stationary regions). There is a constant C and a norm H[F]{0) depending only 
on F and its derivatives on the hyper-surface {0} x (2Af, cw) x 5*^ such that if F is a solution to the Maxwell 
equations (P)-©, then for all t > 0, r, > 1, {0,(p) £ , 

\F{dt + dr„r-'eA)\ <Cr-^/^u+~'H[Fm, 

\Fidt + dr,Al-2M/r)-\dt-dr,))\ + \nr''eA,r-'eB)e^''\ <Cr-' ( ^^JTTT^) ^[^^O), 

|F((1 - 2M/r)-\dt ~ dr,),r-'eA)\ <Cr-\l + \u-\)-' H[Fm. 

If we restrict to u_ < 0, then 

- 2M/r)-\dt - dr,),r-'eA)\ <Cr-\l + \u^\r'^' H[Fm. 

Under the same hypotheses, then for all t > 0, r^ < —1, (6, (j)) G such that > \, 

\F{dt + dr„r-^eA)\ <Cu+-^H[Fm, 

\F{dt + dr,,{l - 2M/r)-\dt - drj)\ + |F(r-ie^,r-ieB)e^^| <Cu+-' H[Fm, 

|F((1 - 2M/ry\dt - dr,),r-'eA)\ <Cu+-' H[F]{Q). 

This gives a decay rate for all components (since, either r > Cu-(- or u_ > Cuj^ in the far region 

> 1). Outside the outgoing light-cone, where < t < — 1, the decay rates are r~^/^, r^^u-^^^^ , and 
r^^u-^"^/^. These are the same rates as can be obtained in Minkowski space using vector field methods, 
with u_ = i — \x\ in Minkowski space. (Faster decay rates can be obtained using conformal compactification 
and other methods.) Approaching null infinity inside the light-cone, with (1 -|- e)r, > t > r.^ > 0, the decay 
rates are r~^/^u^'^^ , r^^u_~-'^/^, and Thus, the decay rates for the first two components are the 

same as in Minkowski space, but the last component decays more slowly than in Minkowski space. The slow 
decay for this component comes from the slow decay rate of t~^ in stationary regions. 

In addition to the Maxwell equations, one can imagine studying the spin 2 field equations. A spin 2 field 
is a (0, 4) tensor with the following symmetries 

WpayS = - WaPyS (4) 
Wal35y = - Wap-iS (5) 

=0 (6) 

=0, (7) 
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and which satisfies the spin 2 field equations 

VW^5c./3 =0 (8) 

^[.W^5\c.0 =0. (9) 

The symmetries of a spin 2 field are similar to the antisymmetry of a Maxwell field, and the spin 2 field 
equations are similar to the Maxwell equations. If the vacuum Einstein equations are satisfied, then the 
Ricci curvature vanishes, and the Weyl curvature satisfies the spin 2 field equations. In K^+^, the spin 2 field 
equations are a good model for the linearization of Einstein's equation about the Minkowski solution, but 
this is not true for the linearization about other solutions. In Cartesian coordinates on the Christoffel 

symbols and the curvature are zero. If one introduces a perturbed metric on Minkowski space and treats the 
Weyl tensor as a tensor field on the original space-time, then the difference between the covariant derivative 
of the Weyl tensor with respect to the perturbed metric and the original metric will be second order in 
the perturbation. Thus, ignoring second order terms, the perturbed Weyl tensor satisfies the spin 2 field 
equations on the original metric. In this sense, the spin 2 field equations are the linearization of the vacuum- 
Einstein equation about Minkowski space. This is the motivation for studying the spin 2 field in [S]. When 
linearizing around a curved space-time, the Christoffel symbols do not vanish, and the linearized Einstein 
equations do not reduce to the spin 2 field equations. More drastically, there is a Buchdahl constraint [26] 
from applying two covariant derivatives, two contractions, and the spin 2 field equations, 

where Riem is the Riemann curvature of the background. On the Schwarzschild manifold, this forces W{T + 
R,T — R,T ± R,r^^eA) to vanish everywhere. This forces the derivative of other components to vanish, 
so that there is only a finite dimensional family of spherically symmetric solution. These are similar to the 
spherically, and non-decaying solutions, which are discussed in appendix \X[ 

Nonetheless, one can ignore the Buchdahl constraint and study the spin 2 field equations. The system 
has been studied formally as a system of transport equations |22j . Here we report that it is possible to use 
a method similar to that we used for the Maxwell equations. It is well-known in the literature that the 
Bel-Robinson can be used to define a conserved, positive-definite quantity from a time-like Killing vector 
in the same way that the energy-momentum tensor can for the Maxwell field. In addition to the conserved 
energy generated this way, one can use a quantity defined in terms of the time-translation symmetry and the 
vector field K. One can again use a method of "spin-reduction" to introduce a "pseudo-Maxwell tensor", 
Fai3 = Wai3-ysT^ , which satisfies the Maxwell equations. One can use control of the pseudo-Maxwell tensor 
to control the trapping terms for the spin 2 field to control integrated norms of the spin 2 field. One can 
then use the symmetries of the Schwarzschild space-time, the field equations, and Sobolev estimates to prove 

pointwise-in-time decay for the non-spherically symmetric components of the spin 2 field. Clearly this is 
pointless, since the dynamics of the spin 2 field are trivial. However, we expect that a similar analysis will 
apply to the genuine, linearized gravity system. The linearized gravity equations are more complicated than 
the spin 2 field equations because there are terms involving the perturbed Christoffel symbols contracted 
against the unperturbed and nonvanishing Weyl tensor. 

In section [21 we introduce several sets of vector fields to provide a simpler notation for discussing the null 
decomposition of the Maxwell field and symmetries. The null decomposition and spinor decomposition are 
essentially equivalent. We estimate energies in section [3] We first review the use of the energy-momentum 
tensor, and then use it to define an energy and a weighted, conformal energy. The growth of the conformal 
energy is controlled by a trapping term which depends only on the zero-weight component. These satisfy 
a simple wave equation, which allows us to bound the conformal energy. In section [4l we use this bound 
and trade Lie derivatives to prove theorem [TJ In section [51 we use the bounds and integration along null 
geodesies to prove theorem [3l In appendix [XI we show that there are no spherically symmetric components 
of the Maxwell field which decay sufficiently rapidly at infinity. In appendix [51 we analyse the wave equation 
governing the zero- weight component using a simplified version of the arguments in [4]. This simplified 
version does not require a spherical harmonic decomposition. 
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2 Notation 



The main purpose of this section is to collect various vector fields and components of the Maxwell field, so 
that the reader can compare the notation used in different places in this paper and elsewhere. 

We begin with some simple notation. We sometimes use vectors as indices on tensors to denote the 
corresponding component. We use the notation g to denote the metric, uj ~ {8, (p) G for the angular 
coordinate, fJ for the angular derivative, and Q, to denote the metric on S'^. Except in the statement of 
results, we use (t, r„ , 9, </>) coordinates, unless otherwise specified. We use X_ to denote the one form generated 
by lowering a vector field X with the metric. We use the insertion operation \x which takes the (0, m) tensor 
A to the (0, m — 1) tensor \xA defined by 1x^(^1, . . . , Ym^i) = A{X, Yi, . . . , Ym^i). We use C to denote an 
arbitrary constant which may change from line to line in a calculation. 

We will prove estimates for t > 0. Because the Schwarzschild manifold is time symmetric, similar 
estimates hold for t S K. In particular, estimates in stationary regions will remain the same, and in ingoing 
and outgoing directions, t and dt must be replaced by —t and —dt- 

2.1 Coordinates, Bases, and Field Components 

Recall the definition of the coordinates r^, 0, 0, m_, and from the introduction. 
We frequently use the coordinate vector fields 

T^du R^dr,, e^de, $=^0, 

and the corresponding normalized vector fields 

From the definition of r*, the definition of R given here is the same as the one given in the introduction. 

We also use null tetrads. In the standard presentation of a null tetrad, the tangent space is complexified. 
A basis {/ox, '^ox, Jti, ™} in which /ox and are (real) null vectors, fh is the complex conjugate of m, 
g(Zox,'T-Gx) = —2, g{m,rh) — 2, and all other inner products between the basis vectors are zero. If X and Y 
are unit vectors orthogonal to /ox and tIox, a null tetrad can be defined by taking m = X + iY. Because of 
this, we will also call a basis {lcx,nc^, X,Y} a null tetrad if g(?cx,?^cx) = —2, g{X,X) = g{Y,Y) = 1, and 
all other inner products are zero. We will generally ignore the distinction between the two definitions of a 
null tetrad. 

One advantage of null tetrads is that they assign weights to certain quantities. Rescaling /ox and riox 
by A and respectively preserves the null tetrad structure, as does rescaling m by e**' (and fh by the 
conjugate, e"*'*). If, under such a change of basis, a quantity transforms as a power of A or of e~"* then the 
corresponding powers are the conformal and spin weights of the quantity. 

We use several null tetrads. The first is the "stationary" tetrad: 

/ 

n 
m 

fh 

Price uses a basis which is parallelly transported along outgoing, radial, null geodesies, j{s) = (s, s + 
C,6o,(j)o) (in the {t,r^,6,(l)) coordinates). The null vectors are 

l^{l-2M/r)-\T + R), 
h=T~ R, 



=f + R={l- 2M/r)-^''^{dt + dr,) 

^f-R={l- 2M/r)-^''^{dt - dr,) 

1 i 
=e + i$ = -de + —^d^ 
r r sm 6 

=6 - i$ = -de "—d^. 

r r sm 6 
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and the angular basis vectors remain the same. To prove theorem [31 in section [SJ we use a basis adapted 
to ingoing, radial, null geodesies. Certain expressions are simplified by using the following coordinate-like 
vector fields. 

L ^dt + dr,=T + R, 
N ^dt - dr, = T-R, 

M=e + ^$ ^dg + ^d^. 

sm a sm 

Christodoulou and Klainerman [5] avoid complexifying the tangent space by using an orthonormal basis 
tangent to the sphere at each point. We use ca and es to denote an orthonormal basis on S^. Thus, r~^eA 
and r~^eB are unit vectors in the Schwarzschild manifold. The indices A,B,... are used for directions 
tangent to the sphere. In summary, we have three null tetrads and a coordinate null basis, 

{Z, n,m,m}, (10) 
{Z,n,m,m}, (11) 

{l,n,r~^eA,r~^eB}, (12) 
{i,7V,M,M}. (13) 

The bases can be used to define the corresponding components of the Maxwell field. The electric and 
magnetic decomposition was already explained in the introduction. Wc now introduce a null decomposition 
and spinor components. These are very closely related but differ in the notation and slightly in the definition. 
The null decomposition consists of two scalars, p and cr, and two 1-forms tangent to spheres, a and a. The 
spinor components arc three complex-valued functions. These are defined in terms of the tetrad in ()12p and 
in Ho]) by 

a{eA) =F{i,eA) (l>i =F{l,m) 

p=^F{i,n) <j)o^]^{F{ln)+iF{fh,m)) 

a^^F{eA,eB)e^^ 
a{eA) =F{n,eA), 0_i =F(r7,, to). 

The spin component index in refers to both the conformal and spin weight. These components are related 

by 

(pi —a{m) 



-1 



=a{m) I'/'-iP =|c^P- 



The null decomposition, a, p, cr, and a, more accurately represents the geometric behavior of the components. 
The spinor notation reveals the spin and conformal weight more easily, simplifies several expressions, and 
suggests connections between the spin wave equation, the spin 1 Maxwell equation, and the spin 2 equations. 
We typically write expressions in terms of the spinor components but think in terms of the null decomposition. 

The spinor components in [22] arc slightly different from the ones we use. Since has conformal weight 
i, replacing the null tetrad (|10p by ([TT]) will take (pi to (1 — 2M /ry/'^fjii. These are the components initially 
used in WB. 
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The spinor components are related to the electric and magnetic components by 



01 + B^) + i{E^ - Sq) 



Certain calculations are simplified by using the null basis (|T3l). We define Maxwell field components 
associated to this null basis by 

$1 ^F{L, M) =r(l - 2il//r)i/20i 

$0 N){1 - 2A//r)-V2 + F{M, M)) =r^(j)o 

=F{N, M) =r(l - 2il//r)i/20_^_ 

To discuss the maximally extended Schwarzschild solution in a neighborhood of the bifurcation sphere, 
it is typical to introduce coordinates 

In the exterior region, these range over J7+ € (0, oo) and U- S (— oo, 0). The coordinates (J7+, U-., 9, 4>) can be 
used in a neighborhood of the bifurcation sphere, and the bifurcation sphere corresponds to (J7+, C/_) = (0, 0). 
The surface t — corresponds to U+U- = — 1. This surface extends through the bifurcation sphere to the 
surface t = in the other exterior region. Since 

e 4M =Cre2AJ (1 — 2M/r), 
on the initial surface t = 0, the coordinate vector fields are 

^ Cri/2eOT(i_2M/r)i/2(T + i?), 

Cr^/^e^ (1 - 2M/rf'^{T - R). 



dU+ 

d 
duZ 



Thus, on the initial data surface and near the bifurcation sphere, the coordinate vector fields and 

-^j— are related to I and h by bounded, nonvanishing functions. If {g^, ^-re used to define a 

tetrad, the corresponding components of the Maxwell field are equivalent to (pi. Since these vector fields 
are coordinate vector fields, they commute. To restrict attention to the region near the bifurcation sphere, 
we will often apply smooth, cut-off functions x<o(m+) and x<o(u-) which are smooth, identically zero for 
u+ > 1 and u- > \ respectively, and identically one for u+ > and u_ > respectively. The vector fields 
X<o("+)73i-'X<o(w-)75n-,0,$ stiU commute. 



2.2 Norms and Lie derivatives 

With the goal of applying derivatives to the components of the Maxwell tensor, we introduce several collec- 
tions of vector fields. Since the vector fields O and $ are not smooth, we use the three rotations of about 
the coordinate axes, 0;. We treat these as vector fields on the Schwarzschild manifold. The collections of 



10 



vector fields we will use are 

T={T, e,}, 
x={i?,T, e,}, 

±={R, f,r-ie,}, 

~ d d 

^={x<o{u+)-g^,X<o{u~)-g^,&i}. 

Since the Schwarzschild manifold is static and spherically symmetric, T generates symmetries of the space- 
time. The normalized vectors in X are used to define the norms of the electric and magnetic components 
of the Maxwell tensor. (The three B^'s can be used to define corresponding components of the electric and 
magnetic field. Taken together these give the angular components.) On the initial data hypersurface t = 

and near the bifurcation sphere (ie, where < 0), the coordinate vectors in X can also be used to define 
the Maxwell field components. 

We now recall some convenient notation for discussing collections of vectors and scalar functions from 
[5]. For two sets of vector fields, A.;, the covariant and Lie derivatives are 

Ca,A2 ={^X^X2\X,eA,}, 

VA1A2 ^{\7x,X2\X, e AJ. 

For two such sets and a (0, 2) tensor A, the components of A with respect to the vector fields are the collection 
of scalar functions 

A(Ai,A2) ={AiXi,X2)\X, e A,}. 

Similarly, for a set of vectors A and a collection of scalar functions {/;}, the derivatives are defined as 

Ca{M - Va{/.} = A{/,} ={Xf\X G A, / e {/,}}. 

This definition holds since the Lie, covariant, and directional derivatives are the same operation when applied 
to scalar functions. For tensor fields, a similar notation can be used to generate collections of tensor fields 
and to consider their components. For example, 

(/:aiA)(A2,A3) ={{Cx,A)iX2,X3)\X, e A,}. 

The same can be defined for iterated Lie or covariant derivatives. 

The norm of a 1-form or a (0, 2) tensor with respect to a set of vector fields is 

\AU^ ^ \AiX,Y)\. 

X,YGA 

The n-derivative norms of a (0, m) tensor with respect to components in Ai and derivatives in the A2 
directions are defined to be 

n 
n 

= E E \{Cx,...Cx,A){Y,,...M\'. 

k=0 Xi,...,XkeA2,Yi,...,Y^eAi 

We note that this notation can be applied equally well with S'^ tangent 1-forms, such as a and a as any 
other forms. 
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3 Control of Energies 

3.1 Summary of the Lagrangian method 

We briefly outline the Lagrangian theory for a general field on a manifold, M . It is assumed that there is 
a scalar Lagrangian i[.T, 0, V0] which is used to define the action, 



S = ( i[.T,0,V(/>]d^a;. 

J M 



If (/) is a minimizcr (or, more generally, a critical point) of the action, then (j) will satisfy the Euler-Lagrange 
equation 

One can then define the energy-momentum tensor from this 

2 r-^^ 



which, by the Euler-Lagrange equation, satisfies 

V"T„^ =0. (14) 
For any vector field X , the generalized momentum vector (^)p and deformation 2-tcnsor ^^^n are 

(^V(r, z) =(VyX)(z) + (VzX)(y) ^""K^p =yc.Xp + v^x„ 

which are related by Stokes' theorem 

an 2 Jq 

This is particularly useful for a Killing vector field, for which ("^^tt = 0. 

For any vector field X, we will define the corresponding energy to be the hyper-surface integral of the 
generalized momentum 



Js 



This depends on the Maxwell field F through the energy-momentum tensor. Frequently, we will be interested 
in t = const hyper-surfaces, for which we define 

Ex[F]it) = / (^)p„di^" = / (^)p„f "(1 - 2M/r)h^dr,d^uj. 

When the deformation tensor vanishes, by integrating over a space-time slab, one gets a conserved quantity: 

Ex[F]{t2) - Ex[F]{h) = 0. 

In applying this integration by parts, we require decay as — > oo, but merely smoothness as r» —> — oo, 
since in the maximal extension of the Schwarzschild manifold, for all values of t, the limit r 2M, tends 
towards the same limiting sphere, the bifurcation sphere. 
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3.2 Quantitative effect of trapping 

The energy-momentum tensor for the MaxweU field is 

Tcp =F^,Fp^ - ^g^pF-'^F^s (15) 
=l{F^^Fp^ + {*F)^,{*F)p''). (16) 

It satisfies ((Ti| and is trace-free 

So^pT"^ = 0. 

FormaUy. one may assume that the MaxweU field is generated from a vector potential Aa hy F — dA, 
Fafj = VqA^ - V/jAa and take the Lagrangian to be L = (1/2)F^'^F^5 = 2{\7^As){V^ A^), in which case, 
the Lagrangian theory for the field A gives the Maxwell equations ((T|) as the Euler-Lagrange equation^ and 
as the energy- momentum, which satisfies . Unfortunately, not all Maxwell fields can be represented 
in this way as an exterior derivative^. However, by direct computation from the Maxwell equation, it follows 
that the energy-momentum tensor in p6|) satisfies p4)) so that Stokes' theorem can still be applied. 

The energy-momentum tensor is strictly positive when evaluated on time-like vectors. We will mainly 
be interested in time-like vectors with no angular components. Since any time-like vector with no angular 
component is a linear combination of I and n, to show the positivity of the stress-energy tensor, it is sufficient 
to compute the components in these null directions. These components are 

T(M")=|0i|2, (17) 

T(U)=|0o|', (18) 
T(n,n)=|0_i|2. (19) 

The Schwarzschild manifold is static, so there is a conserved energy. The energy associated to the 
generator of t-translation, dt, is strictly positive. 



ET[F]{t) =i / (\E\^ + \B\^) (1 - 2M/rydr,d^Lo 

(lap + 2|pp + 2|a|2 + \a\^) (1 - 2AI/rydr^d^Lj 
(101 P + 2|0oP + 10-1 P) (1 - 2M/rydr,d^L0. 



{t}xKxS2 
{tjxRxS^ 



1 

4 
1 

4 J{t}xKxS2 

The corresponding deformation tensor is 

VT = - (1 - 2M/r)-^T® {-^)R+ {I - 2M /r)-^ R® {-^)T, 
=0. 

(We will need to compute deformation tensors later, but, in this case, we could simply have argued that the 
deformation tensor must vanish since dt generates a symmetry of the metric.) From the vanishing of the 
deformation tensor, we have a conservation law 

ET[F]{t) =Et[F]{0). 

This immediately gives an upper bound on the average value of the components of the Maxwell tensor in 
any region bounded away from the event horizon. 

^The other equation, l(2}, holds because = 0. 

*The "magnetically charged solution", F = qg sm{9)d9 A d<f> is not an exterior derivative. See appendix 1X1 
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By applying Lie derivatives, we can get additional conservation laws. If X generates a symmetry and F 
solves the Maxwell equations, then CxF will also be a solution of the Maxwell equations. For each set of 
symmetries and integer k, we have the conserved quantities 

ErlC^om) =£;t[/:^F](0), 
ET[C!^F]{t) =ET[C!iF]{0). 

We now improve these estimates and reveal the effect of trapping, by considering the conformal energy. 
Following earlier work [SI [SJ |U [S] , we let 

K ={t^ +rj^)dt + 2tr^dr, 
^\^{u+^L + uJn). 

We will call this the conformal vector field, but it is also one of the vector fields known as the Morawetz 
vector field. It is an analogue of a vector field used in to prove decay for the wave equation, the 

Maxwell equation, and the spin 2 field. The analogue in R^+" generates a positive quantity, so it is not 
surprising that the same holds on the Schwarzschild manifold. We define the conformal energy to be 

EK[F]{t) = / (^)p„dz." 

J{t}xRx52 




-(1/4) / ((t + + 2{e + r,2)(|p|2 + + (t _ r,)2|a|2) (l _ 2M lr)r^dr,cPuj 

J{t}xKxS2 

= (1/4) / {u+^\4>i\^ + {u+^ + U-^)\4>o\^ +uJ\(t)-i\''){l-2M/r)r^dr,d'u. 

J{i}xRxS2 

In the null decomposition or spinor representation, all the terms in the integrand are non-negative, and, 
inside the light-cone |r*| < (1 — e)i, the coefficients on the Maxwell field components grow like . Thus, 
once wc show that the conformal energy is bounded, there will be decay for the localized field components. 

The following lemma gives an almost conservation law for the conformal energy It states that, to bound 
the conformal energy, it is sufficient to prove sufficiently strong decay in a particularly region bounded away 
from the event horizon. There are two important observations to make from this lemma and its proof: 
(i) an estimate for the two field components Ep^ and Bj^ will control all the field components through the 
conformal charge, and (ii) it is sufficient to control these field components only in a region near the photon 
sphere r = 3M. 

Lemma 4 (Trapping lemma). There is a positive function Xtrap supported in a bounded range of r^ values 
such that if F is a solution to the Maxwell equations (H])-©, then 

EK[F]{t2) - EK[F]it,) < I txtrap + B|) (1 - 2M / r)r' drj'udt 

<[ txtrap\(poWl~2AI/rydr^d^Ljdt. (20) 

J[ti,t2]xRxS2 
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Proof. The deformation tensor for K is given by 



VR^r-'^ {I - 2M/r)g - r-\l - 3M/r){~f® f + R), 
VK + VT + 2tr^yR 

- (1 - 2M/r)-^T ® 2tT + (1 - 2M/r)-^R ® 2r,T 

- (1 - 2M/r)-^T ® 2r^R + (1 - 2M/r)-^R ® 2tR, 
(^V =2trJ"K + U{-f ®T + R®R) 

=4t— (1 - 2A//r)g + 4t f 1 - — (1 -SM/r)) {-f®f + R®R). 
r \ r / 

Because the MaxweU energy-momentum tensor is trace-free, the contraction of the first term against T is 
zero at each point. The importance of r = 3M, where the orbiting geodesies are located, is immediate from 
the second term. The contraction against the energy-momentum tensor is 



(^V„;3T"'' =4t (^(1 - 3M/r) - l) (T^^ - T^^), 
=4i (^(1 - 3M/r) - T(f, h). 
From this, we have the following almost-conservation law 

EK[F]{t2) - EK[F]{h) =2 [ t(l- -(1 - 3M/r)) (^| + S|) (1 - 2M /rydr^Sujdt 

t (^1 - —(1 - 3M/r)^ l0oP(l - 2M/r)r^dr^d^ujdt. (21) 



=2 



[ti,t2]xRxS2 



We refer to 1 — ^(1 — 3M/r) as the trapping term. 

For r 2M, —oo and 1 — 3A//r —1/2, so 1 — (1 — 3Af/r)r,/r is negative. The explicit expression 

for in terms of r is 

fr - 2M\ 

r, =r + 2M log ( j - 3M + 2M log 2. 

Because of the logarithmic term, as r — > oo, (1 — ^(1 — 3Af/r)) = {r — r^,)/r + 0{l/r) < — 2Mlog(r') -|-0(l/r) 
which is negative for sufficiently large r. Since the trapping term has negative limit at ±oo, it is positive 
only in a compact interval. 

We now introduce a smooth, compactly supported function Xtrap which dominates the trapping term. 
This function is chosen to satisfy 4 ((1 — 2,M/r)r^/r — 1) < Xtrap- This gives the desired result. □ 



3.3 Spin reduction 

In this section, we obtain a decay result for the zero-weight component. From the previous section, we 
know this is enough to control energies involving all the components. It is known that the evolution of the 
zero- weight component can be determined from a wave equation without referring to the other components. 
Thus, we can reduce the problem from the Maxwell equations to a wave equation. Since physicists refer to 
wave equations as spin equations and the system of the Maxwell equations as a spin 1 system, we use "spin 
reduction" to refer to this reduction. 

The Maxwell equations can be written as a fairly simple system in terms of the null coordinate bases 
and the corresponding components. This is a result due to Price [22], although, he uses a null tetrad, which 
makes the corresponding expressions look significantly different. By direct computation and application of 
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the Maxwell equations, 



7V$i =M$o(l - 2A//r)r"^ (22) 

L$o =^$1 + cot 6i$i, (23) 

iV$o = -^/*-i -cot6l$_i, (24) 

L$_i = - A7$o(l - 2A//r)r"2^ (25) 

We refer to these as the "Price equations" . 

The cotangent terms appear to be singular; however, if $1 is treated as spherical l-forms, then the 
combination of the angular derivative and the cotangent term can be written simply as 

M$i + cot e**! =(1 - 2Mlr)-^l'^ (diva + icurla) , (26) 
+cot6i$_i =(1 -2M/r)"i/2,^diva-icurla), (27) 

where div and curl arc the spherical divergence and curl. If we'd defined a coordinate based null decompo- 
sition A{eA) = F(L, ca) = (1 - 27\//r)"^/^Q;(eA), then we'd have exactly M$i + cot Q^x = divA + zcurlA, 
and similarly for the other components. One important consequence of this is that the right-hand sides of 
(Pg)) and (gl) are controlled by 



|M$i -t-cot6l$i| + +cot6i$_i| <(1 - 2M/r)-i/2r|F|x 



i,c 



Another important consequence of the Price equations is that the zero weight term satisfies a wave 
equation. From jSH) and ([SS]) . 



LiV$o =(M + cot 6l)Af $0(1 - 2il//r)r"^ 
~d^% = - dl^o + r-2(i _ 2M/r)(-As2)$o. 

If there were an additional (2A//r'^)(l — 2A//r)$o term on the right, then $0 would be a solution to the 
wave equation on the Schwarzschild manifold, V"VQ(r^^$o)- Even in the absence of this term, the previous 
analysis of wave equations is sufficiently general to apply to a wave equation of this form [3]. In fact, the wave 
equation under consideration is simpler than the true wave equatioiH V"Va(r-i u) = 0, and, in appendix 
iBl we provide a stream-lined version of the method from [3] . 

For solutions to a wave equation, there are estimates on the weighted space-time integral we need to 
control for the conformal estimate. If u is a solution to 

~dfu^^df.^u + VL{-As2)u (28) 

with 

VL^r-^l~2M/r), 
then the energy and conformal charge are defined in terms of their densities by 
e=\dtu\^ + \dr,u\^ + Vl\ Wu\^, 
ec =l\it + r)idt + dr,)u\^ + r){dt - drM' + ^(*' + ^^\^ + 

2 J{t}xRxS2 

Ec[u]{t) =- / ecdr^d^uj. 

^ J{t}xRxS2 

'■"If V" VQ(r-in) = 0, then u satisfies -d^u = -d^^u + Vu + VL{-Ag2)u with V = 2Mr-^{l - 2M/r), thus, the equation 
governing <I>o is closer to V'V a{r~^u) = than to V^Vatt = 0. 
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The energy and conformal energy are generated from T and K. (In fact, the conformal energy is generated 
hy K -\-T io provide better control at time t near 0.) Since T generates a symmetry, the energy is conserved. 
As with the conformal energy for the Maxwell field, the conformal energy is not conserved, and the energy 
density near the photon sphere r = 3M must be controlled. The important results for this discussion are 
that, at any time t > 0, fc G N, and any compactly supported function x, 

E[u\(t) =E[um, 
E[JI^u]{t) =E[f^u]{fd), 

Ec[u{t)] <Ec[u{Q)] + CE[^l.u{Q)], 



dr^d^ujdt <E[u]{0), (29) 



[0,oo)xRxS2 (1 + r^'^y 

tx\ JJu\^dr^d^ujdt <CEc[u]{Q) + E[lS.l2u]{Q). 



'[0,oo)xRxS2 

In appendix \^ we exclude spherically symmetric harmonics, so, from dropping the angular derivatives on 
the left-hand side of the previous estimate, 

/ txt,M^{l-'2M/r)r'^dr^d^i^dt<C{Ec[u]{fd) + CE[/\l^u]{fd)). (30) 

J[0,oo)xRxS2 

We apply this with u = $o in the following lemma. 
Lemma 5. If F is a solution to the Maxwell equations and n > 0, then 

(n+l n+5 \ 

k=0 k=0 I 

[ txtrap\<PoWl - 2M/r)r^dr,d^ujdt <C Ek[C^F]{0) + V St[/:^F](0) ) . (32) 

J[0,^)xRxS= Vfc=0 fcO / 

Furthermore, if the normal to S has uniformly bounded below I and n components, then 

I (|(/)i|2 + 2|</)o|2 + |0_i|')(l-2M/r)r2dr,d2tj <C m&^{u^-^ ,u+-'') [ ^ £;k[4-F^](0) + ^ -Bt [4^1(0) ) • 
•'^ \fc=0 fc=0 / 

Proof. Taking $o = u, the energy associated to the Maxwell field F and that of the scalar wave u are closely 
related. From the Price equations and the geometric interpretation of the M + cot 6 terms in (|27p . 

E{u]{t)=ET[CoF]{t), (33) 

E[Al,u]{t)=ET[ClF]{t), (34) 

Ec[u]{t) ^EK[CoF]{t). (35) 

Estimate (l30l) can be written as 



tXtrap|<^o|'(l-2A//r)r2dr,d2^ <C{EK[CoFm + y^ET[C%F]{0)). 

[0.oo)xRxS2 

From this estimate and the trapping estimate (|20p . the estimate (PT|) follows. If one of the derivatives in pi]) 
is in the angular direction instead of the time direction, then it would not be necessary to drop the angular 
derivative in (|30p . and only n and n + A derivatives would be needed on the K and T energies respectively. 

Since the integral of the trapping term has been controlled in the entire of the exterior of the Schwarzschild 
manifold; we have a uniform bound on the integral of on any hyper-surface. If the hypcr-surfacc has 
a normal with uniformly bounded below I and n components, then the integral will be bounded below by 
C(it+^|(/)ip + (it„^ +w+^)|^o|^ + u_^|0_i|^). This provides the final estimate. □ 
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We remark that from equations ([HS]) , we could have bounded the energies Et [^qP] (0 and Ek [CoF] (t) 
by immediately appealing to results for the wave equation. However, we would still need to present the en- 
ergy and conformal energy for the spin wave equation and for the spin 1 Maxwell equation, relate them, 
and present the Price equations. While it would have been possible for us to omit the Lagrangian theory 
and the trapping lemma, this would have removed the motivation for considering the energy and conformal 
energy. 

4 Pointwise decay in stationary regions 

In this section, our goal is to prove L°° decay in regions where 2M < ri < r < r2- We refer to these as 
stationary regions since the range of the radial coordinate does not change in t. Restricting attention to a 
stationary region, the integrand in the conformal energy behaves like times the Maxwell field components 
squared. Since the conformal energy is bounded, the field components decay in L^^^ like t^^. 

Control on radial derivatives is the main thing that we need to improve this from decay in mean to 
pointwise decay. Sobolev estimates can be used to convert L^^^, decay for derivatives into decay. For 
this, we need decay on the spatial derivatives of the Maxwell field. From spherical symmetry, the Lie 
derivative of the Maxwell field in the direction of an angular derivative, Cq.F, also satisfies the Maxwell 
equations and has the same type of decay in mean as F. Since R does not generate a symmetry, the Lie 
derivative in that direction will not solve the Maxwell equations. 

To control the radial derivatives, we use the structure of the Maxwell equations. Using the staticity of 
the Schwarzschild manifold, we can control t derivatives, CtF, in L^^^. In a fixed, compact range of r values, 
the covariant derivatives of the coordinate basis are controlled by finite multiples of the coordinate bases 
again. We are working in where we already control all the components. Thus, we control the difference 
between components of the covariant derivative in a direction and the covariant derivative of the components 
of the MaxwcU tensor (ie, V^F^^X^^rT ~ Va(^)3^X^rT)). 

The notation in subsection l2 . 21 can be used to define a "big-O" notation to estimate the difference between 
two functions depending on position and a tensor field. We say a function of position and a (0,2)-tensor 
field is equal to another such function up to norm terms and in an interval, if, on any bounded interval of 
r* values, there is a constant such that, for any (0, 2) tensor, the difference between the two functions is 
bounded by a constant times the norm of the tensor 



if, for any bounded interval in r*, there is a constant C such that each |/| is bounded by C times the sum 
of the absolute values of the h's plus 0(|A|x) terms. We make similar definitions involving 0{\A\x,k,T)- 

The big-O notation used here is local to compact intervals, which allows us to ignore the difference 
between normalized and unnormalized vector fields. 



This notation allows us to prove theorem [TJ As outlined in the beginning of this section, our strategy in 
the proof is to use the Maxwell equations to trade derivatives along the generators of symmetries for spatial 
derivatives and then to apply the Sobolev estimate. In doing this, we use the big-O notation to estimate 
error terms generated by converting between Lie and covariant derivatives. This allows us to improve our 
decay estimates from decay in mean to pointwise decay. Here, we use the null decomposition and explicitly 
state the norms. Clearly the same result holds for \E\ + \B\ or |a| -I- \p\ + \a\ + \a\. 




\A\^ =0(|A|x), 

l^lx-o(|Ay. 
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Theorem 6. Let 2AI < ri < r2 < oo. There is a constant C^n.r-,) such that if F is a solution of the Maxwell 
equations (U)-©, then for all i G M, r G [?'i,f2]j O'T^d {0,(j)) £ S'^ , 

/ 4 8 \ 

\k=0 k=0 I 

Proof. Using the big-0 notation, wc can control the difference between the derivative of a component of the 
Maxwell field and the corresponding component of the Lie derivative. Since the Lie derivative of any vector 
field with respect to any other is a linear combination of the coordinate vector fields with smooth coefficients, 

Cx{F{X, X)) =iCxF){X, X) + 0(\F\x). 

This process can be iterated, so that 

4(i^(X,X))=0(|F|x,fe,T). 

Similarly for covariant derivatives, 

X(F(X, X)) -(VxF)(X, X) + 0(\F\x). 

We note that if we had applied two symmetry-generating derivatives before making the estimate we would 
have 

CtCt{X{F{X,X))) ^CrCr{{VjiF){X,X)) + 0(|F|x,2,t), 

and similarly with the Lie derivative replacing the covariant derivative in X. 

To control the radial derivative of components which have no R arguments, we use (|2|). 

CRCfCT{F(T,T)) ^CjCjCr{F{T,T)) 
=£f£fWR{F{T,T)) 
=Ct£t{{VrF){T, T)) + 0(|F|x,2,t) 
=CtCt{{VtF){R,T)) + 0(|i^|x,2,T) 
=£tCt\/t{F{R,T)) + 0(|F|x,2,t) 
=£t£t£t{F{R,T)) + 0(|F|x,2,t) 

= Oi\F\x,3,T). 

Similarly, to gain control of component with one radial argument, we apply p]), 

£R£f£f{F{R, T)) <£^£^{{V f^F){R, T)) + 0(|F|x,2,t) 
<£t>Ct((VtF)(T, T)) + 0(|i^|x,2,T) 
=0(|F|x,3,t). 

Since F is antisymmetric, there is no need to control components with two R arguments. 

Control of triple derivative terms, of the form £\£j£j{F{\,'%.)), is sufficient to prove L°° decay. From 
the boundedness of the conformal charge, for any interval [r'i,r2] in the exterior, there is a constant C, such 
that (with and r^^,■2 the values of corresponding to r = ri and r = r2) 

[ t^lFgdr^d^uj <CEk[F]. 

J{t}x('-.1,'-.2)XS2 
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A local, inhomogeneous, 1-dimcnsional Sobolcv estimate 



gives 



{i}x{r.}xS2 



\CTCTiFiX,X))\^dLO <C f 




{t}x(r.i,r.2)xS2 



\Cr£tCt{F{X,X))\^ + \CT£T{F{lL,X))\^dr,d^iu 



\F\'x,3,Tdr*d'^uj 



{t}x(r.i,r.2)xS2 
3 



Now applying a spherical Sobolcv estimate, we have 



F{X,X){t,r^,e,<l))\ <Ct-^ 



J2EK[diF]{t) 



k=0 



3 



2 



By lemma [51 the conformal energy at any time is bounded by the initial data (with extra derivatives). 
This gives, in any stationary region away from the event horizon, that the components decay like t^^. □ 

5 Decay outside stationary regions 

In this section, we prove decay for the field components outside of stationary regions. In Minkowski space 
M^"*""^, it is typical to obtain decay estimates in the regions |a;| < (1 — e)t and > (1 — e)t. Because boosts 
are not symmetries of the Schwarzschild solution, decay in a stationary region is different from decay along 
outgoing curves r» ~ (1 — e)t. Similarly, since there is no reflection symmetry r» — > — , the decay rates in 
the regions r^, > and r» < are different. Thus, we obtain decay in stationary regions, outgoing regions, 
and ingoing regions. 

Most of the decay estimates in this section are proven by considering the energy on ingoing or outgoing 
null hyper-surfaces. We will use and S+_ to refer to ingoing and outgoing null hyper-surfaces on which 
w+ and u_ are constant respectively. We will restrict these to the future t > 0. To make estimates on E+_ , 
we can introduce a parameter ti and an approximate surface which extends along the hyper-surface t = 
from the bifurcation sphere (r* — oo) to the intersection of < = with E+_ , extends along E+_ until t = ti, 
and finally continues onto space-like infinity along t = ti. Since the deformation tensor for T is zero, the 
surface integral of the generalized momentum C^^P along this surface will be the same as the integral along 
t = 0. Similarly, since estimate (|32[) says that the integral over the entire exterior region of the positive part 
of the K deformation tensor is bounded, the surface integral of the generalized momentum (^)p over the 
approximating surface is bounded by the initial data. Dropping the positive contribution from integrating 
along t ~ Q and t ~ ti and taking the limit as ti — > oo, 



A similar argument can be made for . Since and E"^ are null surfaces, we cannot apply the last 
part of lemma [5l 

These integrals can be expanded in terms of the Maxwell field components. 
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To obtain estimates for derivatives tangential to this surface, we can convert to the coordinate based com- 
ponents and apply angular derivatives and the Price equations ([22|) - ([25|) to get 



» 4 8 

/ _ (M+2|iV$i|2(l-2M/r)-V2+u_2|iV$ondu-d'w <C{J2 Ek[CI^F]{0) + J2 M^jF^), (36) 

•^^"+ A:=0 k=0 

. 4 8 

/ ^ (w_2|L$_i|2(l - 2M/r)-V2 + M+2|L$o|2) du+(fuj <C(^ Ek[C^F]{Q) + ^ £;t[4^](0)). (37) 

fc=0 fc=0 

We now prove decay in outgoing regions. 

Lemma 7 (Decay for r» > 1). There is a constant C such that if F is a solution of the Maxwell equations 
(HI)-©, then for all t>0,r^>l, {0, (f) £ , 



4 8 " 



1/2 



\cj^,\<Cr-^l^\u+\-'\Y,EK[C\Fm + Y^ET[C\Fm+ sup ^(r 
\U <Cr-^ f^TT-T^)'^' [Y^EK\ll\Fm^Y.W\m)+ sup Y.^r^l^.f\ , 

/ 4 8 \ 

\<\>-x\<Cr-\\-r\u4)-^\Y^EK\C\Fm^Y.W\Fm+ sup ^(r^/^^oM ■ 
For i < r*, 

/ 4 8 \ 

|0-i|<Cr-i(l + |ii_|)-3/2K^iJ^[4F](O) + ^£;T[4^^](O)+ sup ^(r^/^00' • 

fc^ {0}xR+xS=Y / 

Proof. At any point in the far region, r» > 1, we will integrate along a radial, null ray to prove decay. The 
bounds on the conformal charge give decay for integrals along the null rays. The final end point will either 
be at t = 0, where we already have decay, or at = 0, where we have decay by assumption. In this way, 
each component of the Maxwell field will be written as the sum of two terms, both of which decay. The 
typical null rays which we use are illustrated in figure [2l 

There are a number of simplifications in the outgoing region. We can ignore factors of (1 — 2M/r), since 
the ratio between 1 and (1 — 2M/r) is bounded above and bounded below by a strictly positive number. 
Since > 1, we can ignore ratios of r^/r. There is the ordering > > u_. On outgoing null rays, on 
which u_ is constant, the change in is twice the change in r^, and similarly, on ingoing radial, null rays 
the change in u_ is twice the change in r^. 

The simplest application of our method is for the zero- weight component. First, we prove an estimate 
inside the light-cone, for t > r*. We use a radial, null geodesic from (t, r, , 6, (f>) to (t + r», 0, 9, (p) parametrized 
by U-, 

\^o{t,r,,e,cP)\< [ \N'i>o\du_ + \Mt + r.,0,9,(f>)\ 

1 

— 2j„, 1 I / „, 2| ATif, |2, 

'ci 



1 1 ^ 



The end point decays at a rate of t'~^ evaluated at t' = t + r, . We now integrate over the angular variables 
too (and apply Cauchy-Schwartz, so that the integral in duj is inside the square root). The integral in the first 



21 



(ti +r,i,0) 



{ti - ni,0) 




(0,t2 + r*2) 



Figure 2: Null rays in the outer region, > 0. The angular variables have been suppressed. The null rays 
go from a point either to the initial hypersurface t = or to the stationary region r, = 0. 



term is bounded by the conformal charge as given in ([36]) . The second angular derivatives of F will satisfy 
the same estimate, and we can use the second angular derivatives to control the value of the component, 
through a Sobolev estimate. Hence, 



\J tt}xir,\xS^ ,._n 



< 



{t}x{r.}xS2 

4 

Y,EK[Ci^F]{0) + J2ET[C!^Fm 



1/2/4 8 \ 



\A;=0 k=0 / 



\k=0 fc=0 



Outside the light-cone, where i < r*, we integrate over the curve C4, and the end point value is replaced by 
|$o(0, f + r*, 0, (/))|, which decays like r*'"^/^ = {t + r^)~^/'^ (since (j)o decays like r*'"^/^). Thus, we have 



In the region where < 1, instead of using Ek, we could have used Et, which does not have a vanishing 
factor of u_2. Thus, we may replace (m+ — by (u+ — + 

Now, we prove decay for </>! by again integrating along ingoing, radial, null geodesies again. From any 



-1/2 / 4 8 \ 1/2 
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given point, we integrate along C4 to the endpoint where t = 0, 
0,(^)1 < / \N^i\du_ + |$i(0,i + n, 0,0)1 



< (^J^ u+^\N<^i\'^r^du^^ ' u+- + |$i(0, t + r,,e, (P)\. 

The endpoint wiU be bounded by rj^'^l'^ = + r*)"'^/^. We now integrate in the angular variables, 
differentiate in the angular directions, and apply a spherical Sobolev estimate to get 



l<i>i(i,r„ 0,0)1 <Cu+-V-i/M ^ii;x[4^](o) + ^£;T[4^](o) 

\fe=0 fc=0 J 

+ sup ^(r^/^0.)', 



+ x52 ■ 
4 8 



|0l| <C^.+ -V-3/2 ^i?^[4F](0)+^i?T[4^^](0)+ ^(^5/2^^)2 
\k=0 fc=0 {0}xM+xS^ ,^ ; 

Finally, for 0_i, we integrate along outgoing, radial, null rays on which u- is constant. Inside the 
light-cone t > we take the curve C2 from {t, r^,9, 0) to {t — r*, 0, 0, 0). The estimate is 



^ C-2 



<Cu_-^( / u_2|^$_^|2^2^y^)i ^ |$_i(i-r,,O,6',0)|. 

In the stationary region, the decay rate is also [t')^^ — (t — r*)^^, so the decay rate is 

/ 4 8 \ 5 

|<i>_i(i,r., 0,0)1 <C|u_rM ^i?K[4F](0) + ^i?T[4^](0) , 



\k=0 k=0 
/ 4 8 



I0-1I <C\u^r'r-' [Y^EKidiFm + J2ET[CiiF]{0) 



\k=0 k=0 J 

For t < r^, a similar argument can be made by integrating along C3, with the value at the other end point 
being (f>_i(0, r» — t, 6, 0), where we have faster decay, 

|<I>_i(i,r., 0,0)1 <C|u-r3/2 [^£;;,[4F](0) + ^i?T[4^](0)) +\u^\-'/' sup ^(r^/'0O', 

/ 4 8 \ 5 

I0-1I <CK|-'/V-M^ii;;,[4F](0) + ^i?T[4^](0)+ sup ^(r5/20,)2 . 

\t^0 t^O {0}xR+xS=Y / 

Again, in the region < 1, we can use Et instead of Ek to get a better bound when u- vanishes. □ 
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We now turn to proving decay in the "near" region, < 0. Since for any fixed interval 2M < ri < r < r2, 
we can apply the results from section |4l the main purpose of the following lemma is to prove estimates which 
are uniform in r so that they can be extended to the event horizon. Note that extends smoothly to the 
event horizon and is an affine parameter for tangential geodesies. 

Since the vector fields T and K vanish on the bifurcation sphere, the boundedness of the associated 
energy allows rapid divergence of the (normalized) energy density near there. Not surprisingly, this is not 
sufficient to control the Maxwell field. In light of this, it is somewhat surprising that the energies associated 
with T and K are sufficient to prove decay for the correctly normalized components of the Maxwell field 
tensor corresponding to (pi . 

As explained in the introduction, the correctly normalized basis for stating results on or near the event 
horizon is 

dt + dr„ il~2M/r)-\dt-drJ, r-^CA, r-^es. 

We can equally well replace r~^eA and r~^eB by m and rfi or by G) or $. For large this is the "correctly 
normalized" basis, because it is the result of parallelly transporting the original, normalized basis on the 
initial data surface, < = 0, along ingoing null geodesies to reach the event horizon. 

The method used in the previous lemma gives decay rates of and for $i, <I>o, 

and <I>_i respectively. The functions <I>i and $0 are correctly normalised (except for bounded factors) as 
r 2M, but (1 — 2i\//r)^^<l>_i is the correctly normalised component in this region. For $0, we prove 
a different preliminary decay rate and then use the divergence theorem to obtain a rate of For the 

correctly normalised, negative-weight component, we use the decay for $0i ^ transport equation, and 

an integrating factor to get u+~^ decay. We note that the vector field (1 — 2M/r)^^{T — R), which is a 
smoothed version of the vector field y in [D] can be used to prove boundedness for this component without 
using a transport equation. 

Lemma 8 (Decay for < 0.). There is a constant C such that if F is a solution of the Maxwell equations 
then for all i > 0, < 0, (6*, 0) G 5"^ such that u+ > 1, 

\F{dt + dr,,e)\ + \F{dt + dr,,n <cu+-^ (^i?A'[4^^](o) + 5^i?T[4J^](o)) , 

\fe=0 fc=0 / 

\Mt, r.)\ <C\Fidt + 9.., (1 - 2M/r)-\dt - 9.J)| + |F(e, l>)| 

( 4 8 3 \ 5 

^ Ek [CtiF] {0) + Y,Er [4^] (0) + E I^^l (0) , 
A:=0 fc=0 fe=0 / 
/ 4 8 3 \ 

F((l - 2M/r)-\dt ~dr,),Q + i^) <Cu+-^ \ y2EK[C^TFm + Y.^T[diFm + Y,Ef[ClFm 

\fc=o fe=o fc=o / 

Proof. The first part of this proof is similar to that of lemma [71 The main difference is that we must track 
factors of (1 — 2M / r) carefully, but we may ignore factors of r since it is bounded above and below by positive 
constants. When tracking factors of (1 — 2M/r), we use (1 — 2M /r)(j)i{t' ,rj , 9' , tp') to denote the value of 
(1 — 2M/r)(j)i at [t' ,rj 9' ,(j)') even if an unprimed set of coordinates is in use simultaneously. Since we are 
only considering > 1, the ingoing, radial, null rays from any point will hit the stationary region r* = 0. 
This is illustrated in figure [3] 

Integrating along an ingoing, radial, null geodesic, C5, from (t, r*, 0, (j)) to [t + 7%, 0, 0, (j)), in the same way 
as in the proof of lemma [71 we have 



|$i(t,r*,0,0)| < / |iV$i|d?/_ + |$i(t + r„O,0,0)| 

<(/ {l-2M/r)duS][ I \N^i\^u+^r'^{l-2M/r)-^du-] u+-^ + \^i{t + r^,Q,9,(j))\. 
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(t + r„0) 



Figure 3: Null rays in the inner region r* < 0. The angular variables have been suppressed. The curve C5 
goes from a point in the stationary region = to an arbitrary point in the regions t > 0, r* < 0, u+ > 1 
along an ingoing, null, radial geodesies. 



The integral of (1 — 2M/r) with respect to du- is, up to a factor of 2, the same as the integral of (1 — 2M/r) 
with respect to dr*. By a change of variables, this is the integral of 1 with respect to dr. Thus, the 
contribution from the first integral in the first term is bounded by a constant. After integrating in the 
angular variables, the second integral is bounded by the conformal energy of the angular derivatives of F. 
The second term is the value of the component in the stationary region, so it decays like (t + r^)~^. Applying 
the angular derivative and Sobolev estimate argument from the previous lemma. 



I $1 1 = + a,. , e + i<l) I <c?.+ - M ^ [4^] (0) + 51 [4^] (o) • 

\fe=0 k=0 / 

This proves the first result. 

For the zero weight component, we first prove an intermediate result for {1—2M/ r)^/^(/)o. From integrating 
along a surface of constant t, we have 

{l-2M/r)^/'-\u+\^^lit,r,,e,<P)=~ [ dr, ({1 - 2M/ry^^\u+\^<Pl) drj 

J{t}xlr,.,0]x{e}x{4,} ^ ' 
+ (1 - 2il//r)l/2u^2^2(^ 

The integrand can be estimated by dropping negative terms and applying the Cauchy-Schwartz inequality 
as 

-dr, ((1 - 2M/ry/V\M') = - ^(1 - 2M/r)i/2^«+2|</>o|' - 2(1 - 2M/ry^'uM' 

- 2(1 - 2M/r)^/\t+2K((/)oa,..^o) 
<C {u+^l - 2M/r)\cj,o\'' + u+^\drM^) ■ 

Thus, for u+ > 1, the integrand is controlled by the conformal energy using the standard Sobolev estimate 
argument by 



r / 4 8 

/ dr, ((1 - 2M/ry/^\u+f^l) drj <C V i^^i^^KO) + V i?T[4^^] 

J{t}x[r,.0]x{9}x{4.} ^ ^ \^ ^ 



(0) 
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t 




(0,r, -t) (0,n + f) 

Figure 4: The region Q,(^t^r,)- 



The end point can be controlled by the stationary decay result, 

/ 4 8 N 

(1 - 2Af/r)i/2u+2 |ci,o {t + r,,0,d,(f>)\^ <C[Y,Ek [O^F] (0) + J] i?T [OiF] (0) 



^A;=0 fe=0 



SO that 

/ 4 8 \ 

\/c=0 fe=0 / 

This estimate can now be improved. For a given point (i, r*, 0, (^), consider the two-dimensional surface 

17(4 ,,^) ={(i',r,',6',</)) : t' > 0,t' + rj < t + n,t ~ < t' ~ rj <t + i\}. 

This is illustrated in figure [3 Applying the (two-dimensional) divergence theorem with the vector field 
(Af$o)^, we have 

<i>oit,r^,e,(l))-<i>a{t + r ^,0,6,(1))+ N<i>adr^=-2 LN<^Qdr^dt. 

J{O}x[r.-t,-t-r.]x{e}x{0} •'^''■[t.r,) 

From the wave equation for $0, we have 



-2/ r-\l-2M/r){~/^s2)<^fidr^dt. 

We now estimate the terms on the right. The first is bounded by the stationary decay result. The second 
is an integral in the initial data surface i = 0, so it can be controlled by integrals of the initial data. We have 

I / N<j)odr,\<i [ {1^2M/r)^^^drA if {1 - 2M/r)-^/^\N(f>o\^drA . 

J[r,-t.-t-r,] \J[-oo,-t-r,] J \J [-00 ,-t-r,] J 

The first integral is bounded by 

/ (1 - 2Mlrf/^dr, <C(1 - 2M/r)i/2 < c'e(-*-'-*)/2A^ < Cv 

•J \—oo.—t—rA 
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Since (1 — 2M/r) ^^'^N and d/dU- differ only by smooth functions of r, the second can be estimated by 

[ (1 - 2A//r)-i/2|iV0opd'-* <C f Ij^MH^ " 2M/r)'/^dr,. 

Since the stationary tetrad and the one based on the J7_, 6, (j) coordinate system differ only by smooth 
functions of r, 

/ {l~2M/r)-^l^\NM''dr. <C [ |^ f ^^(U) + F(e, \\l - 2M/r)'^'dr, 

J[-oo-t-r,] J[-oo,-t-r,] (^U - \Z ) 

+ ij'^f Jr) + ^<«- - ^'W"*v 

Since coordinate vector fields commute, 

/ (1 - 2M/r)-^/^\N(f,o\^dr^ <C [ V \C''_j^F\l{l ~ 2M/r)^^^dr,. 



k=0 



The same argument could have been applied to the second angular derivatives of F, which could have been 
used in a Sobolev estimate. This would have lead to 

/ (1 - 2M/r)"i/2|iV0opdr, <C V i?^[4F](0). 

Thus, the integral along the initial time slice is bounded by 

3 \ 1/2 



J[r,-t-t-r,] 



Finally, we estimate the integral over ri(t,r.) by breaking it into two parts, ri(t ,,,)^ = ^(t,r.) H {t > 2|r*|} 
and fl(t,r,}g = ^{t,r,} n {< < 2|7%|}. In fl(t,r,), 

r-^{l - 2M/r){-As2)<i>odr^dt\ 
< sup ((l-2M/r)i/4(-A52)$o) / r'^il - 2M/r)^/^dr^dt 

r-^{l-2M/r)dr^dt \ ( / r-2(l - 2Af/r)|A52$opdr,di ) 

On the first line of the right-hand side, the supremum term decays like by the intermediate result, and 

the integral term is uniformly bounded by the exponential decay of (1 — 2M/r) with respect to r^,. In the 
second line, the second integral is bounded by estimate ([29]). and the first integral is bounded by (1 — 2M/r) 
evaluated at the point {t\rj) = ((2/3)u+, (l/3)u+), and hence decays faster than u+~^ . Combining all 
these results gives 

/ 4 8 3 \ 

\k=0 fc=0 fc=0 / 
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We begin our analysis of <I>_i with an intermediate decay result, using the same sort of simple argument 
as was used for $i. Integrating along an outgoing, radial, null ray, we have 



|$_i(i,r,,6',(/))| <M_-i (^y"(l-2A//r)du+^ (^J \L'i>^i\^uJ{l ~ 2M/r)-^du+^ + (t - r,, 0, 61, 0) 

/ 4 8 \ 

\A;=0 fe=0 / 

Since wc arc working in the inner region, r* < 0, there is the estimate w+ < u_, and 

/ 4 8 \ 

\k=0 k=0 ) 

To obtain stronger estimates, we will need to integrate along outgoing null geodesies starting near the 
initial data surface. We work with the [/+ and [/_ coordinates to control the correctly normalised, negative- 
weight component near the bifurcation sphere. If bounded initial data is posed on the surface U+U- = — 1, 
which corresponds to the union of the t = surfaces in the two exterior regions, then, at least in in some 
small neighborhood of the bifurcation sphere, the components of the Maxwell field with respect to the 
smooth {U+,U-,9,4>) coordinate system must remain bounded by a multiple of their initial value. This is 
essentially a Cauchy stability result, as was used for the wave equation [S]. Thus, in some sufficiently small 
neighborhood of the bifurcation sphere, in the exterior, the outgoing component is bounded by 

3 
fc=0 

In particular, we can pick a it+Q <S such that on the hypersurface u+ = u+q where t>0, 

3 

(l-2Af/r)-i|$_i| <cY^Ef[Cp']{Q). (38) 

fc=0 

We now use the decay for <I>o and <I>_i to prove a stronger estimate. From the Price equation (|25p . we 
have 



Since 1/2M — IMjr^ vanishes linearly at r = 2Af, it is bounded by C(l — 2M/r). Integrating along an 
outgoing geodesic starting on u+ = it+Q at (io, ?'*oi 0) ^^id going to (i, r*, 6*, </>), we have 



*0: ' 



The geodesic along which we integrate is illustrated in figure [5l We break the integral into two pieces, with 

one going from u+g to M+/2 and the other going from w+/2 to From the boundedness and decay for $o 
and $-1, the integral is bounded by 

» / 4 8 3 

•^^e \fc=0 A;=0 fe=0 
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Figure 5: An outgoing null ray from the hypersurface uj^ = it_|-Q to the point under consideration. The angular 
variables have been suppressed. The region under consideration contains a portion near the bifurcation sphere 
and is near the event horizon. Therefore, the t and coordinates are not used. Instead the hypersurfaces 
t = 0, r = 2A/, and u+ = u+q are indicated. 



By estimate ((38)) . the end point, (1 — 2M/r) ^$_i(to, r^g, 6*, 0) is bounded. Thus, 

|(l-2M/r)-icI>_i(i,r,,0,0)| <C(e-«+/«^^ + ^.+ -l) ( ^ £^[4^1(0) + ^ i?T[4^](0) + ^ [^FKO) ) 

\fc=0 k=0 k=a / 

/ 4 8 3 \ 

+ Ce(-"++"+o)/4M ^ (0) + [C\F] (0) + 5] i?^ [ClF] (0) 



X 

Vfc=0 fc=0 k=0 



/ 4 8 3 \ 

\A:=0 A;=0 fc=0 / 



This provides the desired decay. 



Combining the two lemmas in this section, we have theorem [31 



□ 



Remark 9. Finally, we provide a simpler sufficient condition for the initial data to have finite norm. 

The energies Et[F] and Ek[F] are the integrals of the field components squared, (jff, against the weight 
(1 — 2M/r)r'^ and rj^{l — 2M/r)r'^ respectively. Thus, if the field components are hounded everywhere and 
decay at least as fast as r~'^/2+'^) (for e > as r ^ oo, these energies will be finite. Thus, if the Maxwell 
field and its first eight derivatives decay at this rate, then 

4 8 

j2EK[ctiFm+j2ET[c!iFm 

k=0 k=0 

will be bounded, and the result of theorem [3| will apply. 

Similarly, if the field components decay at least as fast as r~'-'^/2+'), then sup|o}xM+xS2 
he trivially bounded, so that the results of lemma^will give decay in the far region > 1. 

The energy Ej, is the integral of the field components squared, , against the measure {l — 2M/ry^^r^dr^. 
Thus, if the components are bounded as — > — oo, then this energy will he finite. Since (for < 0), the 
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vector fields in X are coordinate vector fields extending in a neighbourhood of the bifurcation sphere, they are 
smooth vector fields. Thus, if F and its first three derivatives with respect to this collection of smooth vector 
fields have finite components, then 

fc=0 

will be finite, and the results of lemma will apply. 

This verifies the footnote to theorem [7] that boundedness and r-^s/^+e decay for F and its first eight 
derivatives on the initial data surface t = is sufficient to prove the decay results in this paper. The same 
remark applies to the result in theorem\^ 

A Exclusion of the non-radiatable mode of the Maxwell field 

In this section, we show that if the MaxweU field has finite conformal charge, it has no sphericaUy symmetric 
part and explain why this is a physically reasonable assumption. The absence of dynamic, spherically 
symmetric components is well-known in the literature [22]. In K-'^+'^, since the electric and magnetic fields 
are divergence free, the spherically symmetric component of the Maxwell field is always zero. While there 
are spherically symmetric solutions on the Schwarzschild manifold, we show that they have no dynamics, 
since these solutions are constant in t. There is a two-parameter family of such solutions described by the 
central electric and magnetic charge. These solutions do not vanish on the event horizon and decay like l/r^ 
at infinity, so they fail to be in the finite conformal energy class we consider. 
The Maxwell field F can be written as 

F(t, 61, 0) =p(t, 61, - 2M/r)At A dr* 

+ r{l~2M/rf/^dt^uJo{t,r^,6,(j)) 
+ r{l~ 2M/rf/^dr^ A Wi {t, r* , 6*, </)) 
+ <T{t,r^,e,(t)yVL, 
*F{t,r^,9,(l)) =-cr(I-2Af/r)dtAdn 

- r(l - 2M/r)i/2dt A (^s^ Vi 
-I- r(l - 2M/r)i/2dr* A (>i=s2)cjo 
-t- pr'^n, 

with p and a scalar functions, with {*s^) the Hodge dual on S^, with loq and uji in f7^(S'^) for each value 
of t and r*, and with f2 the standard volume form on S^. We have used a stationary, instead of null, 
decomposition, so ujq and uJi appear instead of a and a. The scalars p and a arc the standard ones from the 
null decomposition. 

We first remind the reader that there is no cj G 0^(5^) with ds^cu = {*g2)Co. If there were one, we could 
write this condition in coordinates: 

ll> —ujgdO + uj^d4>, 
(^9,4, - =Co sin(6'). 

Let f{0) — J^^gi (^cj>d(f). Since uj is smooth, / is continuous on [0, 27r] and vanishing at and tt (since the 
integral is over a single point in in these cases). From duu ~ {*s^)Co, we have f'{0) = Cosin(6'), and 
/(■"') ~ /(O) = /o C'o sin(6')d6' > 0. Thus, the condition that / vanishes at and tt implies Cq = 0. 
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The Maxwell equations are (taking the components orthogonal to various 1-forms) 



Equation Orthogonal 1-form 



dF 



-At : 



-At : 



dr, : 



dr, : 



















t(1 
t(1 
t(1 
t(1 



2M/ry/^dt A ds2UJo + {dt{(Tr^))dt A n 
iM/rf^dr^ Ad52(*52)wo + (<9^. (/9r^))dr, ^n 
2M/rY''^dt ^ds2{*s^)^l + {dt{pr'^))dt ^^l 



Since dt and dr* are spherically symmetric, these can be projected onto the I = spherical harmonic 
(equivalently, we can contract on the 5^ volume). Since there is no / = component for the 1-forms wq, wi, 
(*52)a;o and 32)001, we find, 



Thus, there is no dynamics is the I = mode, since the t derivative is always zero. These solutions do not 
decay sufficiently rapidly to have finite conformal energy. 

The exclusion of these spherically symmetric solutions is physically reasonable. Physically, the solutions 
represent a perturbation of the Schwarzschild black hole to a charged Reissner-Nordstrom solution, not 
an external perturbation by radiation. Price refers to these spherically symmetric solutions as the "non- 
radiatable" modes, since the solutions in this two parameter family are static. Since the Maxwell equations 
are linear and commute with angular derivatives, the spherically symmetric component does not couple to 
the other components, so it will not affect the rest of our analysis to eliminate the spherically symmetric 
components. In analogy with the theory of solitons, we might think of the Reissner-Nordstrom solutions as 
a manifold in the space of solutions to the Maxwell-Einstein system. In this case, the decoupled Maxwell 
equations with I > correspond to linearized perturbations from this manifold, whereas perturbations with 
I = correspond to linearized perturbations along the manifold of stationary solutions. 

B Analysis of the wave equation 

We now prove decay estimates for solutions to the wave equation ((28)) . For this equation, there is also an 
energy and conformal energy, which we use in our analysis. As with the Maxwell field, we must control the 
trapping of u near the photon sphere to control the growth of the conformal charge. We do this with a local 
decay estimate and employ light-cone localization to obtain a local decay estimate of the full strength we 
require. The arguments and results of this section are a slight modification of the argument in [J, only, in 
this case, the situation is simpler. 

For the wave equation on the Schwarzschild manifold, previous analysis [2 |4l |9] has required a decom- 
position onto spherical harmonics. On each spherical harmonic, the wave equation can be treated as a 
one-dimensional wave equation with an effective potential. The main estimate uses a vector field, 7, which 
points away from the maximum of the effective potential. In the case of the geometrically defined wave equa- 
tion, the location of these maxima depend on the spherical harmonic parameter, and 7 has been modified to 



5..(ai=or') =0 
dtipi^or'') =0. 



From which it follows that the I = components are given by integration constants qe and qs, 
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fit each spherical harmonic. The equation (|28p is simpler, with the maxima of the effective potential always 
at =0. Thus, a very minor modification of the previous analysis allows us to make the estimate without 
using a spherical harmonic decomposition. 

Since the potential Vl ~ -fji^ ^ 2M/r) is real- valued, we may analyse the real and complex parts, which 
each satisfy separately. Thus, we may assume our solutions are real- valued. 

We use the method of multipliers to analyse Although it maybe possible to introduce a Lagrangian 
formulation and an energy-momentum tensor, we do not do so because this is not the geometrically defined 
wave equation V"Vct(r~"'^M) so the Lagrangian for this system is quite artificial, because there would be a 
confusion between the energy-momentum for u and that of the full Maxwell field, and because some of the 
energies would require correction terms. The essence of the method is to choose a "multiplier" (a differential 
operator), apply it to the function u, multiply by the equation, and integrate by parts. The most useful 
differential operators are typically those given by the vector fields from the Lagrangian method. 

We begin by recalling the energy and conformal energy, which were defined in section [S] Conservation of 
energy follows simply from the method of multipliers with the multiplier T = dt- Multiplying (|28p by dtu, 
we find 

=(dtu)(-dfu + u+ 4(1 - 2M/r)As2u) 
= -ldt(,\dtu\^ + \dr,u\'' + VL\ Wu\^) + drAdtudr,u)+ y ■ {VLdtu Wu). (39) 
Integrating over a space-time slab gives conservation of energy: 

E[u]{t) - E[u]{0) =0. 
Using the method of multipliers with K gives 

Ec[u]{t2) - Ec[u]{ti) = / 2t{2VL + r^V[)\ fu\^ dr^d?ujdt 

J[ti,t2]xExS'2 

< / iXtrapl Wufdr^d^ujdt. (40) 

J[ti,t2]xRxS2 

This is similar to the estimate for the Maxwell equations. For this analysis, it is useful to introduce an 
energy localized inside the light cone. We let 

Emin = / edr^d^Lu. 

J{t}x[-(3/4)t,(3/4)t]xS2 

We need a variety of Hardy estimates. 

Lemma 10. Ift> 1, xh is a non-negative function which is positive in some open set |r, | < t, and a > 0, 
then if f :Rx S"^ is a smooth function, and u : [ii, x K x -S*^ — > K is smooth with t G [ti, t2\ and u{t) ~ f , 
then 



dr^d'^u <CEn^in[u]it), 



{t}x[-(3/4)t,(3/4)t]x52 (1-f 7-,2) 

;dr,d'u;<C f ,f''^^C +XH\ffdr,d'LO 



{t}x[-(l/2)t,(l/2)t]xS2 (l + |r'*|)"+2 7{t}x[-(l/2)t,(l/2)t]xS2 (1 + f*) 
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Proof. We start working with a > 0. When > 0, 



(l + r,)"+i ' '(l + ?v)"+2 



(l + ?%)«+2 '-(a + l)2y(, + a + 1 

Since, for any exponent /3 > 0, (1 + r*)'' is equivalent to (1 + rj^)^/"^ on [0, cxd), the powers of (1 + r*) can 
be replaced by (1 + r*2)i/2. By symmetry, the same result holds on (— r^xi 0]. Since {\ + rJ^)~^ is uniformly 
equivalent to (1 + (r* — r*o)^)~'^ foi' ''*o in a finite interval, the |/(0)|2 term can be replaced by |/(7'*o)P in 
any fixed interval. 

For a = 0, we take r^^ = (3/4)i. By integrating the estimate over r^g with r^Q in (f/2,3/4), where 
is strictly positive, and then integrating over the angular variables, we find 

/ TrT^dr,d^u:<CE^M{i)- 

J{t}x[-(3/4)t,(3/4)t]xS2 l-LH-r* j 

Similarly, taking r<,i — (l/2)t and a > 0, for any non-negative function, xh, which is positive in some open 
set inside |r* | < t, 

J{t}x[-(l/2)t,(l/2)t]xS2 U + nj ^ J{t}x[-(l/2)t,(l/2)t]xS2 U-|-|r*|j 

□ 

We now prove a local decay estimate to control the trapping terms. To do this, we use a radial multiplier 
7 in terms of a weight 

7 = .95r. + {dr,g)/2. 

Assuming that the weight g is a function of the t and r» variables only, we have, 

—dt{2u^u) = — dt{2ugdr,u + u{dr,g)u) 

= - dr, {9{dr,uf + {dr,u){dr,g)u - 14g()i7u)2 + - f . {VL{Wu)gu) 

+ 2{dr,g){dr,uf - _ {d,^VL)g\ fu? 
— 2ug{dr„u) — u{dr,,g)u. 

We use the notation 

Ej[u]{t) — / u{'-fu)dr„d^Lij. 

J{t}xRxS2 

Integrating over a space-time slab, 

-2E^[u]{t2) + 2E^[u]{ti) = / {2{dr,g){dr,u)^ - idr,VL)g\ JJu\^ 

J[ti,t2]xRxS2 ^ 

—2ug(dr^u) — u{dr,g)u)drt,d^ Lodt. (41) 
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Taking X[-i/A.?,/A] to be smooth, non-negative, compactly supported in [—3/4,3/4], and identically 1 on 
[—1/2, 1/2], b to be a sufficiently small parameter to be chosen later, and a e (1, 2], we set 

XLC =X[-3/4,3/4](y), 

We now expand and group the terms on the right of (|41]) . 

2txLcidr J){dr,uf - iXLC (5^. Vl)^! fu\^dr,Sujdt (42) 



-2E^[u]{t2) + 2E^[u]{ti) = / 



t2]xRxS2 



2 



^Xlc ' '"l" — dr^Suidt (43) 

[ti,t2lxRxS2 ^ 

2tg{dr,Xhc){dr,ufdr^d'^ujdt (44) 

[ti,t2]xRxS2 

i(3(a,..XLc)(9'.5) + 3(c'2^XLc)(9r.ff) + (a.^.XLc)^)^^^^^^^^ 

[ti,t2lxRxS2 ^ 

(45) 

2ug{dr,u) + u{dr,g)udr^d^ijjdt. (46) 

[ti,t2]xRxS2 

The terms on the right in lines (|44|) - (P5|) can be estimated by the local energy. (In these calculations, 
remember that t^^ < r^,^^ and inverse powers of t arise from differentiating xlc-) 



/ 2tg{dr, 



Xi.c){.dr,ufdr^d'^uj 



<C I {dr,u) dr^d LO < CEn^in, 

J{t}x[-(3/4)t,(3/4)t]xS2 

/ \mdr,xi.c){dlj) + i{dl,xi.c){dr.g) + {dlxi.c)~g)^\dr,d^u 

J{t}xRxS2 ^ 



}xRxS; 



<C [ ^ , ] .M '^dnd^Lu < CE^i^, 

J{t}x[-(3/4)t,(3/4)t]xS2 t + r*| 

\2ug{dr,u) + u{dr,g)u\drifd^ijj 

1 

+ \c)r,u\^ + - — ;■ 

'{t}x[-(3/4)t,(3/4)t]xS2 



<C I \u\'^ + \drM^ + --^\u\^dr,d^uj 



<Ci?niin- 

The left-hand side can be estimated similarly by 



E^{t) <tC [ 
J{t} 



|m|2 + \dr,u\'^ + - — ]—^\u\'^dr^d^uj 

x[-(3/4)t,(3/4)t]xS2 1 + |r. 



<CtE„ 



The two terms on the right appearing in line (|42p are clearly positive, since g is increasing and was chosen 
to go from negative to positive at the same value of as — V^- To control the term in by the terms in 
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, we note that 

dl^~g =a{cT + 1)52(1 + - 



{t}x[-(3/4)t,(3/4)t]xS2 ^ J{t}x[-(3/4)t,(3/4)t]xS- 

We divide the range of integration into two pieces and use the estimate t/l?', | < C when > t/2. From 
this, 



{t}xRxS2 ^ 2 J{t}x[-(l/2)t,(l/2)t]xS' 



2 J{i}x{|r.|e[(l/2)t,(3/4)t]}xS2 



Applying the Hardy estimates, we find 



/ 

J{t} 



- txLC-^\u\'^d.r^(fuj 

J{t}x[-(l/2)t,(l/2)t]x52 



2 

Taking 6 suSiciently smah, we can dominate the integrand by half the terms in line (|42p . Thus 

1 
2 



tXLC (2(9rJ)(ar.u)2 - t{dr,VL)g\ ^up) di\(fujdt < - 2£;^|*J + C / E^in{t)dt. 

[ti,t2]xKxS2 -/ti 



The same estimate holds for Cq.u, so that summing over the components, we have 

\ I txLC {2{drJ){dr, ](7u)^ - ixLC (5,-. Vl )g | u| 2) dr^d^iodt < - 2E^[](/u]\ll +C T E^-^,,[J/u]{t)dt. 

^ J[ti,t2]xRxS2 Jti 

The left-hand side controls the trapping term by the Hardy estimate. Since the trapping term controls the 
growth of the conformal charge, 

Ec[u]{h) <£;cM(0) + I - 2E^[fuf^-\ + C f ' E^,^[Wu]{t)dt, 

<CEc[u]{0)+ sup {tE^in[fu]it)) + C E,^in[](7u]{t)dt. 
te[o,t2] Jo 

Since is conserved, there is an immediate linear bound on the conformal charge. By applying the 

Caucliy-Schwartz estimate and integration by parts both twice, we can make the estimate 



EnuniW^m <E[A's2u]{t) 



This allows us to make a self-improving estimate. From the linear bound, the conformal energy can't grow 
faster than i^/**, and the t^/* implies a uniform bound. Thus, 

Ec M (t) <CiEc [u] (0) + EiAl^u] (0)), 

txl y/ufdr^d'^LO <C{Ec[u]{0) + E[Al2u]{0)). 

[0,oo)xRxS2 
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Applying the same argument with the factors of t and xlc dropped, so that g{t,r^) = g{r^), we find 



Acknowledgement 

This project was started when the author visited UC San Diego. The author would like to thank J. 
Sterbenz for his hospitality and for contributing several valuables ideas. 

References 

[1] A. Bachelot. Operateur de diffraction pour le systeme de Maxwell en metrique de Schwarzschild. C. R. 
Acad. Set. Paris Ser. I Math., 312(1):93~96, 1991. 

[2] P. Blue and A. Soffer. Phase space analysis on some black hole manifolds. (math.AP/0511281). 

[3] P. Blue and A. Soffer. The wave equation on the Schwarzschild metric. II. Local decay for the spin-2 
Regge- Wheeler equation. J. Math. Phys., 46(1):012502, 9, 2005. 

[4] P. Blue and J. Sterbenz. Uniform decay of local energy and the semi-linear wave equation on 
Schwarzschild space. Comm. Math. Phys., 268(2):481-504, 2006. 

[5] D. Christodoulou and S. Klainerman. Asymptotic properties of linear field equations in Minkowski 
space. Comm. Pure Appl. Math., 43(2):137-199, 1990. 

[6] D. Christodoulou and S. Klainerman. The global nonlinear stability of the Minkowski space, volume 41 
of Princeton Mathematical Series. Princeton University Press, Princeton, NJ, 1993. 

[7] M. Dafermos. The interior of charged black holes and the problem of uniqueness in general relativity. 
Comm. Pure Appl. Math., 58(4):445~504, 2005. 

[8] M. Dafermos and I. Rodnianski. A note on energy currents and decay for the wave equation on a 
schwarzschild background. (arXiv:0710.0171vl). 

[9] M. Dafermos and I. Rodnianski. The red-shift effect and radiation decay on black hole spacetimes. 
(gr-qc/0512119). 

[10] M. Dafermos and I. Rodnianski. A proof of Price's law for the collapse of a self-gravitating scalar field. 
Invent. Math., 162(2):381-457, 2005. 

[11] J. Dimock and B. S. Kay. Classical and quantum scattering theory for linear scalar fields on the 
Schwarzschild metric. II. J. Math. Phys., 27(10):2520-2525, 1986. 

[12] F. Finster, N. Kamran, J. SmoUer, and S.-T. Yau. Decay rates and probability estimates for massive 
Dirac particles in the Kerr-Newman black hole geometry. Comm. Math. Phys., 230(2):201-244, 2002. 

[13] F. Finster, N. Kamran, J. SmoUer, and S.-T. Yau. Decay of solutions of the wave equation in the Kerr 
geometry. Comm. Math. Phys., 264(2):465-503, 2006. 

[14] F. Finster and J. SmoUer. Linear stability of the schwarzschild black hole under electromagnetic and 
gravitational perturbations, (gr-qc/0607046). 

[15] D. Hafncr. Sur la thcorie de la diffusion pour I'cquation de Klein-Gordon dans la metrique de Kerr. 
Dissertationes Math. (Rozprawy Mat.), 421:102, 2003. 

[16] D. Hafner and J. -P. Nicolas. Scattering of massless Dirac fields by a Kerr black hole. Rev. Math. Phys., 
16(1):29-123, 2004. 




36 



[17] S. W. Hawking and G. F. R. Ellis. The large scale structure of space-time. Cambridge University Press, 
London, 1973. Cambridge Monographs on Mathematical Physics, No. 1. 

[18] W. Inglcsc and F. Nicolo. Asymptotic properties of the electromagnetic field in the external 
Schwarzschild spacetime. Ann. Henri Poincare, l(5):895-944, 2000. 

[19] B. S. Kay and R. M. Wald. Linear stability of Schwarzschild under perturbations which are nonvanishing 
on the bifurcation 2-sphere. Classical Quantum Gravity, 4(4):893-898, 1987. 

[20] C. W. Misner, K. S. Thorne, and J. A. Wheeler. Gravitation. W. H. Freeman and Co., San Francisco, 
Calif., 1973. 

[21] C. S. Morawetz and W. A. Strauss. Decay and scattering of solutions of a nonlinear relativistic wave 
equation. Comm. Pure Appl. Math., 25:1-31, 1972. 

[22] R. H. Price. Nonspherical perturbations of relativistic gravitational collapse. IL Integer-spin, zero-rest- 
mass fields. Phys. Rev. D (3), 5:2439-2454, 1972. 

[23] R. H. Price and L. M. Burko. Late time tails from momentarily stationary, compact initial data in 
Schwarzschild spacetimes. Phys. Rev. D (3), 70(8):084039, 6, 2004. 

[24] J. V. Ralston. Approximate eigenfunctions of the Laplacian. J. Differential Geometry, 12(1):87-100, 
1977. 

[25] T. Regge and J. A. Wheeler. Stabihty of a Schwarzschild singularity. Phys. Rev. (2), 108:1063-1069, 
1957. 

[26] J. Stewart. Advanced general relativity. Cambridge Monographs on Mathematical Physics. Cambridge 
University Press, Cambridge, 1990. 

[27] S. A. Teuksolsky. Rotating black holes: Separable wave equations for gravitational and electromagnetic 
perturbations. Phys. Rev. Lett., 29:1114 - 1118, 1972. 

[28] B. F. Whiting. Mode stabihty of the Kerr black hole. J. Math. Phys., 30(6):1301-1305, 1989. 

[29] F. J. Zerilli. Gravitational field of a particle falling in a Schwarzschild geometry analyzed in tensor 
harmonics. Phys. Rev. D (3), 2:2141-2160, 1970. 



37 



